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We provide a quasi-one-dimensional model which can support a pair-density-wave (PDW) state, 
in which the superconducting (SC) order parameter modulates periodically in space, with gap¬ 
less Bogoliubov quasiparticle excitations. The model consists of an array of strongly-interacting 
one-dimensional systems, where the one-dimensional systems are coupled to each other by local 
interactions and tunneling of the electrons and Cooper pairs between them. Within the interchain 
mean-field theory, we find several SC states from the model, including a conventional uniform SC 
state, PDW SC state, and a coexisting phase of the uniform SC and PDW states. In this quasi-ID 
regime we can treat the strong correlation physics essentially exactly using bosonization methods 
and the crossover to the 2D system by means of interchain mean field theory. The resulting critical 
temperatures of the SC phases generically exhibit a power-law scaling with the coupling constants of 
the array, instead of the essential singularity found in weak-coupling BCS-type theories. Electronic 
excitations with an open Fermi surface, which emerge from the electronic Luttinger liquid systems 
below their crossover temperature to the Fermi liquid, are then coupled to the SC order param¬ 
eters via the proximity effect. From the Fermi surface thus coupled to the SC order parameters, 
we calculate the quasiparticle spectrum in detail. We show that the quasiparticle spectrum can be 
fully gapped or nodal in the uniform SC phase and also in the coexisting phase of the uniform SC 
and PDW parameters. In the pure PDW state, the excitation spectrum has a reconstructed Fermi 
surface in the form of Fermi pockets of Bogoliubov quasiparticles. 


I. INTRODUCTION 

The pair-density-wave (PDW) state is a superconduct¬ 
ing (SC) state of matter in which the Cooper pairs have 
a finite momentum. Due to the finite momentum carried 
by the pair, the SC order parameter is modulated peri¬ 
odically in space. The PDW state has recently received 
attention because it can explain the layer decoupling ob¬ 
served in the cuprate La2- x Ba a: Cu04(LBC0), the orig¬ 
inal high T c superconductor, at the x = 1/8 anomaly^ 
At this doping the T c of the three-dimensional material 
is suppressed to temperatures as low as 4K, where the 
Meissner effect first appears and the system is in a three- 
dimensional d-wave SC phase. In contrast, away from 
x = 1/8, the SC T c in LBCO is about 35K. In spite of 
the low T c of the uniform d-wave SC state at x = 1/8, 
in this doping regime the CuO planes are already super¬ 
conducting for a range of temperatures well above T c £& 
With this in mind, Berg and collaborators^ showed that 
this phenomenon can be explained if the CuO planes are 
in an inhomogeneous (‘striped’) SC state, the PDW state, 
in which charge, spin and SC orders are intertwined with 
each other. In this state the SC order parameter os¬ 
cillates along one direction in the CuO planes, and the 
average of the SC order parameter is zero in the CuO 
planes. 

The PDW state is similar to the traditional Larkin- 
Ovchinnikov (LO) stated where the Cooper pairs have 
a non-zero center of mass momentum, which in the LO 
proposal is due to the presence of an external Zeeman 
field, which breaks the time-reversal symmetry explicitly 
(for a recent review of the LO state see Ref. 3). How¬ 
ever, the occurrence of the PDW SC state does not nec¬ 


essarily require to have a system in which time-reversal 
symmetry is explicitly broken, nor it does require time- 
reversal symmetry be spontaneously broken either. Since 
it was proposed as a candidate competing state to the 
uniform d-wave SC order^ the PDW state has been 
studied extensively. A Landau-Ginzburg (LG) theory of 
the PDW state provides a simple explanation of much of 
the observed phenomenology of La2-:rBa x Cu04,— ^ and 
of La2- a; Sr a ;Cu04 in magnetic fields. An outgrowth of 
these phenomenological theories is a statistical mechani¬ 
cal description of the thermal melting of the PDW phase 
by proliferation of topological defects which yielded a rich 
phase diagram which includes, in addition to the PDW 
phase, a novel charge 4e SC state and a CDW phased— 
More recently, Agterberg and Garaudi^ showed that it is 
possible to have a phase in which a uniform SC and PDW 
order parameters coexist in the presence of a magnetic 
field. 

The microscopic underpinnings of the PDW state are 
presently not as well understood as the phenomenolo- 
gies. Nevertheless, it has been shown that this state 
can appear in different regimes of several models. In 
the weak coupling limit, such a state appears naturally 
in two dimensions (2D) inside an electronic spin-triplet 
nematic phased Also at the mean field level, Lee^ 
found that it is possible to have a PDW state in his 
model of ‘Amperian’ pairing^ and that the PDW state 
can explain the pseudo-gap features found in the angle- 
resolved photoemission experiment^ In a series of pa¬ 
pers, Loder and collaborator s 17 ! 18 found that a PDW su¬ 
perconducting state is preferred in a tight-binding model 
with strong attractive interactions (although the critical 
value of the coupling constant above which the PDW is 
stable is presumably outside the range of validity of the 
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weak coupling theory). Similarly, PDW states with bro¬ 
ken time-reversal invariance and parity have been found 
recentl y 19 ' 20 in a ‘hot spot’ model, which also requires a 
critical (and typically not small) value of a coupling con¬ 
stant. On the other hand, in one-dimensional systems 
(ID) the PDW state has been shown to describe the SC 
state of the Kondo-Heisenberg chain,— i 21 ' 22 and a phase 
of an extended Hubbard-Heisenberg model on a two-leg 
ladder— 2 We showed recently that the PDW state ap¬ 
pearing in these two ID models is actually a topological 
SC which supports Majorana zero modes localized at its 
boundaries— - 

There has also been considerable recent effort to de¬ 
termine if the PDW state occurs in simple models of 
strongly correlated systems. Variational Monte Carlo 
simulations of the t — J and t — t' — J model on the 
square lattice at zero magnetic field near doping x = 1/8 
found that the uniform d -wave SC state is slightly fa¬ 
vored over the PDW state— Sr— Corboz and coworkers,— 
using infinite projected entangled pair-states^ (iPEPS), 
found strong evidence in the 2D t — J model that the 
ground state energies of the uniform d -wave state and 
the PDW state are numerically indistinguishable (within 
the error bars) over a broad range of dopings and pa¬ 
rameters. This last result indicates that these strongly 
correlated systems do have a strong tendency to exhibit 
intertwined orders and that the PDW state occurs more 
broadly than was anticipated.— 

In this work, instead of following a conventional weak 
coupling approach to the PDW states, we will take an 
alternative path which has the physics of strong corre¬ 
lations as its starting point. Rather than starting from 
a true 2D system, we will consider a quasi-one dimen¬ 
sional model consisting of weakly coupled (each strongly- 
interacting) ID systems. In the decoupled limit we can 
solve each ID system non-perturbatively using bosoniza- 
tion methods.—— We will follow a dimensional crossover 
approach that has been used with considerable success by 
several authors— — We will consider a generalization of 
the model used by Granath and coworkers^ in which 
there are two types of ID subsystems: a set of doped 
two-leg ladders in the Luther-Emery (LE) liquid regime 
(which has a single gapless charge sector and a gapped 
spin sector) and a set of ID electronic Luttinger liquids 
(eLL) with both a gapless charge sector and a gapless 
spin sector. Although the interactions between LE liq¬ 
uids and between LE and eLL liquids will be treated by 
the interchain mean field theory (MFT) (see, e.g. Carl¬ 
son et al.—), the intra LE and intra eLL interactions 
are treated essentially exactly using bosonization. We 
will make the reasonable assumption that the interac¬ 
tion between the electronic Luttinger liquids leads to a 
crossover to a full 2D (anisotropic) Fermi liquid (see, e.g. 
Ref. [36|). In this fashion, this approach allows to ac¬ 
cess the strong coupling regime of a strongly correlated 
system using controlled approximations. In this approx¬ 
imation the resulting superconducting T c is a power law 
in the interchain coupling and not exponentially small as 


in the usual weak-coupling limit (such as BCS approach). 

The main departure of the system that we consider 
here from previous studies of models on this type is that 
we will allow for the Josephson couplings between the 
LE liquids to have either positive or negative signs. A 
negative sign induces a tt phase shift between two neigh¬ 
boring LE liquids. It was shown by Berg et air that two 
superconductors that are proximity coupled to each other 
through a ID weakly doped Hubbard model have a broad 
regime of parameters (in particular, doping) in which the 
effective Josephson coupling is negative. Here, in order to 
incorporate this physics, we will introduce a set of Ising 
degrees of freedom mediating the interactions between 
the LE liquids which emulate different doping profiles 
of the electronic Luttinger liquids. This feature will al¬ 
low us to consider the interplay between uniform (s-wave 
or d-wave) superconductivity with PDW superconduct¬ 
ing states and coexistence phases, resulting in complex 
phase diagrams. 

We note that inhomogeneous SC states such as the 
PDW are generally accompanied by a subsidiary charge- 
ordered state, a charge-density-wave (CDW). The period 
of the CDW is twice the period of the PDW or equal to 
the period of the PDW depending on whether this is a 
pure PDW state or if it is a state in which it coexists 
with an uniform SC state. The CDW order parameters 
which describe these states are composite of the PDW 
order parameters with or without the uniform SC or¬ 
der parameter. The general occurrence of charge-ordered 
states as subsidiary orders of an inhomogeneous SC state 
has been emphasized by several authors— iL£ i 1 :L i 14 i 19 The 
same should hold in the case of the SC states that we 
study here. 

The experimental consequences of the PDW 
states have been discussed extensively in the recent 
hteraturei^&ii' 14 ->^2r— (including papers by one of 
the authors) and we will not elaborate further on this 
questions here. Instead we will focus on microscopic 
mechanisms behind these inherently strongly interacting 
states. 

The paper is organized as follows. In section [TT] we de¬ 
fine our model and summarize our notation for bosoniza¬ 
tion in ID. In section m we develop the interchain MFT 
and discuss the results for the self consistent equations. 
In section[|V]we study and discuss the quasiparticle spec¬ 
trum of the phases, emergent from this quasi-ID system, 
for the various PDW and uniform SC states found in sec- 
tion llin In section[V]we summarize other possible phases 
that could arise in this model using a qualitative scaling 
dimensional analysis. We finish with our conclusions in 
section eh 


II. THE MODEL 

The quasi-ID model, schematically presented in the In¬ 
troduction, consists of two different types of ID systems. 
One of them is a conventional ID electronic Luttinger 
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liquid (eLL) in which both the spin and charge degrees 
of freedom are gapless. The other type, however, is a ID 
system with a spin gap, i.e., it is a ID Luther-Emery liq¬ 
uid (LE). The presence of the spin gap in the ID system 
will bias the full array of ID systems to a strong tendency 
to a SC state. 


A. ID Systems and Bosonization 


Before we define in detail the quasi-ID model, we start 
with a short summary on those ID liquids and their de¬ 
scription using bosonization. This material is standard 
and can be found in several textbooks, e.g. Ref. (31 j. Here 
we will only give some salient results and set up our def¬ 
initions (and conventions) that we use in later sections. 

We start with a ID eLL which has a gapless charge 
sector and a gapless spin sector. The low-energy Hamil¬ 
tonian is written in terms of the set of the bosonic fields, 
{</> o ,0 o } where a = c,s labels the charge and spin sec¬ 
tors respectively. These fields satisfy canonical equal¬ 
time commutation relations 


[<t>a(x'),d x O b (x)\ = iS a:b S(x' - x ) (2.1) 

The effective Hamiltonian for the eLL is 

K a (d x e a ) 2 +, 

a 

( 2 - 2 ) 

in which K a (again with a = c, s) are the Luttinger pa¬ 
rameters for the charge and spin sectors, and v a are the 
characteristic speeds for the charge and spin excitations 
of the liquid. The parameters K c , K s , v c and v s are deter¬ 
mined by the microscopic details of the model. However, 
for a system with spin-rotational invariance, the resulting 
SU(2) symmetry restricts the value of the the spin Lut¬ 
tinger parameter to be K s = 1. In this continuum and 
low-energy limit, we can decompose the electronic field 
operator in terms of two slowly varying components, with 
wave vectors near the two Fermi points ztkp 

4=tM*) Ra(xy kFX + L„(x)e~ ikFX , (2.3) 
yja 


H e LL[(t>a,0a] = Y. Y 


where a is the ultraviolet cut-off (typically the lattice 
spacing), and where a = ± denotes the spin of the elec¬ 
tron. Here the fermionic fields R a {x) and L a (x) are the 
right- and left-moving components of the electron field 
Vv, which are slowly varying in space relative to the 
Fermi momentum Uf- 

The right- and left-moving fields can be written in 
terms of the bosonic charge fields (f> c and 9 C , and the 
spin fields (j) s and 9 S , as follows 


_ F<J iy/Tr/2(9r.+oO,+<l>r.+ir<l>*) 

V2™ 

_ Fa / 2(8 c+aSs—d c.—a<t>a) 

\Zh ra 


(2.4) 


The anticommuting Klein factors, F a , ensure the 
fermionic statistics between the right and left moving 
fermions R a and L a . 

Next we consider a spin-gapped Luttinger liquid, or 
LE liquid. At energy scales below the spin gap A s , the 
spin sector can be ignored. Hence, the low-energy Hamil¬ 
tonian contains only the charge fields <p c and 9 C , and it 
is given by 


Hle[4>c e c ] 


Vc 

2 


K c (d x 9 c ) 2 + ^(d x ^ c ) 2 

c 


(2.5) 


Since the spin sector has been effectively projected-out, 
we will keep only the charge sector of the LE liquid and 
drop its c label. 

In the LE liquid, all interactions represented by op¬ 
erators that are not spin singlets are irrelevant (and, in 
fact, with effective scaling dimension infinite). This fact 
strongly restricts the types of interactions between LE 
systems and eLL systems. In this case the only fermion 
bilinears that need to be considered in the LE liquid are 
the order parameter of the charge-density-wave with mo¬ 
mentum 2k f (CDW) 

OcdwOt) ~ e~ 2%kpx R\[x)L a {x) + h.c., (2.6) 


and the (Cooper) pair field spin singlet superconducting 
order parameter 


A(i) ~ R a (x)(ia v ) a/3 Lp(x) + (R <->■ L ). (2.7) 

Hence the coupling to the LE liquids in the quasi-ID 
model should involve only the two operators listed above. 
We note that the suppression of the spin operator and the 
power-law correlation for the SC operator make the LE 
liquids the natural building blocks for the quasi-ID SC 
state. In contrast, the eLL has other observables that 
need to be considered, including a spin triplet pair field, 
the 2 k F spin-density-wave (SDW) ‘Neel’ order parame¬ 
ter, the right and left moving spin current operators, and, 
in tunneling processes, the electron operators. 


B. Quasi-ID Model 

Given a set of independent ID LE and eLL systems 
that were described above, we now define and discuss the 
full quasi-ID model. The model consists of an array of 
ID systems shown in Fig. [I] Each unit cell of the array 
consists of one LE system, labeled by A, and one eLL 
system, labeled by B. Hence we introduce the bosonic 
fields {<f> n ,A, 9 n ,A } representing the charge fields in the LE 
chain of the n-th unit cell and also {4> n ,B,a, 9 n ,B,aj (with 
a = c, s) representing the charge and spin fields in the 
eLL chain of the n-th unit cell. Furthermore we assume 
that the filling of the type B system (an eLL chain) is 
close to half filling, i.e., k^ ss 7 t/ 2 and Kc^ ~ 1/2. 
Also the spin rotational symmetry in the B systems is 
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FIG. 1. (color online) The array of LE systems and eLL 
systems. The A-type LE systems are represented by the solid 
(red) line. The B-type eLL systems are represented by the 
dashed (blue) line. Each unit cell consists of one A-type 
and one B-type systems. Here the black filled dots repre¬ 
sent electrons, (a) The conventional Josephson coupling Jab 
in Eg. (12.91) (b) The conventional Josephson coupling Jaa in 
Eg. (12.91) (c) Splitting a Cooper pair in A system into the 
neighboring B systems J' AB in Eo. (12.9ll . 


assumed to be unbroken and thus = 1. We further 

assume that the Fermi momenta of the systems A and B 
are incommensurate to each other. 

In the limit in which the LE systems and the eLL sys¬ 
tems are decoupled from each other, the effective Hamil¬ 
tonian of the array is described by the sum of Ea. (12.21) 
and Ea. (12.51) for each system, and has the form 

Hq = ^ (^H LE [4>n,A, 0 n ,A] + H LL [(j) n ^B,a, 0n,B,a])- 

n£Z 

( 2 . 8 ) 

This decoupled limit is an unstable fixed point and the 
system will eventually flow to the quasi-ID or 2D fixed 
points under the introduction of the coupling between 
the ID systems. We will show that the PDW state, as a 
quasi-ID fixed point, will emerge from certain couplings. 

Following the work of Granath et al.fi* we write down 
the possibly relevant local interactions terms between the 
ID systems. They are given by 

H ' = / dx { ~ tnn i,< 7 + h.c.] 

n ^ a 

H - JBB &B ,j ' &B,j +1 
— •JAA,j[^A i j + \^A,j + h.C.] 

— Jab[^ b + h.c.] 

+ JAB[ A A,j(i’B,j,r^Bj-l,l + V'Bd-l.tV’Bdd) + h ' C '. 

(2.9) 

To simplify the analysis, in this paper we will not con¬ 
sider the possible existence of spin-ordered phases (i.e. 


spin stripes or SDWs) although these are clearly seen in 
La 2 - 2 ; Ba a ;Cu 04 which is the material where the PDW 
state is most clearly seen. We are mainly concerned 
about the SC states in which the spins do not play much 
role, and thus we ignore for now the antiferromagnetic in¬ 
teractions in the discussion. The antiferromagnetic cou¬ 
pling between the eLL chains can also be included in a 
relatively straightforward extension of the present work. 
Following Ref. [3(| we have ignored the possible CDW 
couplings between chains. In general the scaling dimen¬ 
sions of the CDW operators become less relevant in the 
presence of forward scattering interactions between the 
chains^ so they can be neglected. If the interchain CDW 
couplings were to become relevant we would have bidirec¬ 
tional charge order. In this paper we are only exploring 
states with unidirectional charge and superconducting or¬ 
der. 

In Eci. (12.9D . the operator A aj(x) ~ 

’ 4 , A,j,a{x)(i<r y ) al 3 il’A,j, 0 (x) represents the density of 
the spin-singlet Cooper pair of the system A, and 
A B,j(x) is that of the system B. The effective coupling 
constants Jab and Jaa, j are the conventional Joseph¬ 
son coupling between the two neighboring A systems, 
representing the hopping process of the Cooper pairs 
(see the (a) and (b) in FigQ]). On the other hand, the 
local term J' AB represents the breaking of a Cooper pair 
in an A system which puts the two single electrons into 
the nearest neighbor B systems (see the (c) in FigQ] for 
a diagram of the process). 

In the Hamiltonian H ', Ea. (12.9D . the most relevant 
term is the electron tunneling term, with coupling 
strength t B b , between two nearest-neighbor B systems. 
Under this perturbation, the decoupled B systems flow 
to the 2D Fermi liquid fixed point, which, in turn, be¬ 
comes coupled to the superconducting state emergent 
from A systems.— Due to this dimensional crossover of 
the B systems it is difficult to apply the conventional 
interchain MFT to analyze Ea. (12.91) . In order to make 
progress, we ignore at first the B systems, as the first 
order approximation to the problem, and perform the 
interchain MFT only with the A systems, which embod¬ 
ies the strong-coupling nature of the superconductivity 
emergent in the quasi-ID models. We should stress that 
in the A systems, there are no electron-like quasiparti¬ 
cles due the existence of the spin gap which leads to a 
fully gapped 2D SC phases when the coupling between 
the A systems are turned on. Technically speaking, we 
solve first for an array of B (eLL) systems coupled by t B B 
and for an array of A systems coupled only by Jaaj in 
Ea. (12.91) . and take J' AB and Jab as perturbations. At 
this level of the approximation, the emergent SC state 
is determined by the Josephson coupling J A Aj and the 
subsequent SC state of the full system follows by proxim¬ 
ity effect between the A and the B subsystems. This was 
the strategy used by Granth et al..— The main difference 
between this work and that of Granath and coworkers is 
the inclusion of an additional, Ising-like, degree of free¬ 
dom in the coupling between the A systems, as we already 
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discussed in the Introduction. 


C. Coupled LE Systems 

It is clear that Jaa will determine the nature of the 
emergent 2D SC state from the quasi-ID model. More 
precisely, the spatial pattern of J a a. j determines that of 
the Cooper pair. For example, if the Josephson coupling 
Jaa.j in Eg. (12.91) is uniform and positive, it is clear that 
the uniform spin-singlet SC will emerge. However, in 
the strongly-correlated quasi-ID system, the Josephson 
couplingjTAAj may not be always uniform and positive. 
In Ref. (4j, the Josephson coupling between two systems 
with an intermediate chain (which is close to the insulator 
phase) has been calculated by a numerical density ma¬ 
trix renormalization group (DMRG) method and it was 
found that it can be negative, i.e., forming a 7r-Josephson 
junction between the two A systems. From this, it is not 
difficult to imagine that there might be more complicated 
patterns, depending on the microscopic details, than the 
uniform 7r-Josephson junction. 

To reflect this physics and to consider a broader pos¬ 
sible patterns of the Josephson coupling Jaa ,j, we in¬ 
troduce a phenomenological Ising degree of freedom (jj 
which can change the magnitude and possibly the sign of 
effective Josephson coupling. This Ising degree of free¬ 
dom can be regarded as a local change in the doping of 
the intervening B system between two neighboring A sys¬ 
tems. In this sense the Ising degree of freedom should be 
regarded as reflecting the tendency to frustrated phase 
separation of a doped strongly correlated systeim 45 ! 46 To 
this effect, we consider the following interaction between 
the Ising degrees of freedom and the LE liquid 

ffint — 0AA ^ ' COs[\/27r(0A,i ^.4,2+1)] 

i 

-J'aaY. O'iO'i -)_i COs[J 2ir{0 A,i 9 a, i- t-l)] 

i 

+ ^isi„gM, (2.10) 

in which we write Jaa,i in Eg. ( 12 . 91 ) as Jaa,i — Jaa - 
J'AA a i a i+ 1- The factor ~ cos[v / 27r(0i — 0»+i)] is the 
Josephson coupling between the LE systems because of 
A(x) ~ 

In Eci. (l2.10l) . the Ising interaction Hamiltonian 
Rising[o»] is assumed to have several phases depending 
on the parameters in Rising and temperature, e.g. para¬ 
magnetic phase (<Ti) = 0, and various symmetry-broken 
phases. In this paper, we further assume that the Ising 
variable a j orders at a much higher temperature (or en¬ 
ergy scale) than the spin gap A s in the LE liquid. Hence, 
we ignore any correction to the Ising variable due to the 
fluctuations of the SC states emergent from LE liquids. 
To simplify the analysis we have assumed that the Ising 
variables are constant along the direction of the ID sys¬ 
tems and are classical (i.e. we did not include a trans¬ 
verse field term). The first assumption is not a prob¬ 


lem since we will do mean field theory assuming that 
the resulting modulation (if any) is unidirectional. More 
microscopically we will need to assume that the Ising 
model has frustrated nearest and next nearest neighbor 
interactions along one direction only. This is the so- 
called anisotropic next-nearest-neighbor Ising (ANNNI) 
model which is well known to have a host of modulated 
phases4£ Similar physics, with a rich structure of pe¬ 
riodic and quasi-periodic states, is obtained from the 
Coulomb-frustrated phase separation mechanism , 46 1 48 
In what follows we will not specify the form of Rising 
and assume that its ground state is encoded in a specific 
pattern of order for the Ising variables. In this picture 
an inhomogeneous charge-ordered state occurs first (and 
hence has a higher critical temperature) and this pattern 
causes the effective Josephson couplings to have an “an¬ 
tiferromagnetic” sign (i.e. n junctions) Nevertheless, 
as noted in Ref. 3, once the PDW state sets in there is 
always a (subdominant) CDW order state with twice the 
ordering wave-vector as that of the PDW. 

The symmetry breaking patterns that we to study are: 
i) Uniform configuration (a) = ±, Vi, ii) Staggered con¬ 
figuration ( <n) = (—1)% and iii) Period 4 configurations 
(which will become clear soon below HTC 41) . Thus, when 
the Ising variables order and spontaneously break the 
translational symmetry, the effective Josephson coupling 
between the different A systems will be modulated too. 
For concreteness, throughout this work, we w r ill take Jaa 
and J'aa be positive. This condition is not necessary 
and the following arguments can be easily extended to 
the other signs of Jaa and J'aa- 

We will start by analyzing the ground state (or the 
mean field (MF) state) of the LE systems coupled to Ising 
variables. We will do this for different configurations of 
the Ising variables and see what are the possible phases 
that arise in the system of coupled LE liquids. 


1. Ising Paramagnetic Configuration 

Before proceeding to the symmetry-broken phases of 
the Ising variable, we first briefly comment on the case 
with the paramagnetic phase of the Ising variable cq. 
In the Ising paramagnetic phase, we first note that 
J^cqcq +1 cos[V2n(9i — #i+i)] is effectively zero at the 
level of mean field theory and can be ignored. Thus 
Eq. ( 12.10D will become at the low energy 

—Jaa ^ ' cos[\/27r(0A,i — 9a,i+ i)] + • • ■ (2-11) 

i 

in which ■ • • are the terms generated by integrating the 
fluctuations of the Ising variables in the paramagnetic 
phase, e.g., ~ cos[2-\/27r(0A,i — 9a, i+ i)]> which is strictly 
less relevant than — Jaa cos[\/27r(#A,2 — 9a, i+ 0] appear¬ 
ing in "Hint- It is well-known that Eo. (l2.11D induces an 
uniform 2D superconducting statei 34 ' 40 
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2. Uniform Ising Configuration 

We now analyze the simplest case with (a/) J 0, where 
all the <jj have the same value, cq = a = ±. In this case 
lit is just given by 

Hint = — + J'aa) COs[\/2tt(9a ,i — 0A,i+l)] 

.-J T ^cos[v^(^-^ + i)] (2.12) 

i 

The system of coupled LE systems can be treated in in¬ 
terchain MFT, where all the systems are in phase, since 
Jt > 0. In this case we just have a uniform SC state 
in the direction perpendicular to the systems A j = A, 
where A includes the spin gap and the MFT value for 
(cos V27 t6a,%)- We will show in the following section how 
to compute the value (cos VZkOa,/)- Thus, this is the 
same phase as in the Ising paramagnetic case but with a 
larger value of the effective Josephson coupling. 


3. Staggered (Period 2) Ising Configuration 

Let us now consider cq = (—1)L In this case Hint is 
given by 

Hint = -{Jaa ~ Jaa) cos[\/27r(6U,i - 0a,%+ 1)] 

i 

.-^E c °s[^(^-^ + i)] (2.13) 

i 

Again, the system of coupled LE systems can be treated 
in interchain MFT. However, we need to be careful about 
the sign of 5J. If SJ > 0 the SC order parameter in 
all the systems are in phase. It is important to em¬ 
phasize that although all the systems are in phase as 
in the uniform Ising configuration, the expectation value 
(cosv27 T0A,i) is different in both cases, since as we will 
see below, it explicitly depends on the coupling between 
the systems, in this case Jt or 5J. On the other hand, 
if SJ < 0 the phase of SC order parameter has a n phase 
shift between nearest neighbors. In the former case we 
just have a uniform superconducting state in the direc¬ 
tion perpendicular to the systems, while in the second 
case we have a PDW state A aj ~ (—1)L There is a 
direct transition from the uniform SC state to the PDW 
SC state at Jaa/ Jaa = 1- In this simple period 2 Ising 
configuration there is no room for coexistence between 
the uniform SC and the PDW state. 


4. Longer Period Ising Configurations 

We can generalize the phases obtained with period 2 
Ising configurations to cases with longer periods of the 
Ising variables. For instance for a period 4 of the Ising 


variables, 

••• ,t,t,, ( 2 . 14 ) 

the effective Josephson couplings will have a period 2 
modulation. In this case we will find either an uniform 
SC state or a period 4 PDW SC state, but no coexistence 
phase. 

However, we will see that for Ising configurations with 
period n, with n > 2, we can have a richer phase diagram, 
including a coexistence phase if n > 3. For example, for 
a period 3 structure of the Ising variables, the allowed SC 
state is a coexistence phase, whereas for period 8 with the 
following spatial pattern of the Ising degrees of freedom 

(2-15) 

we will find either a coexistence phase with period 4 or 
a PDW SC with period 8. It is straightforward to gen¬ 
eralize this to more intricate configurations of the Ising 
variables. 


III. INTERCHAIN MFT ON THE LE SYSTEMS 

Keeping the quasi-ID model of the previous section in 
mind, we now solve the coupled LE system problem us¬ 
ing the interchain MFT. In this section, we generalize the 
works of Lukyanov and Zamolodchikov,— and Carr and 
Tsvelib^ to the patterns of the Josephson coupling be¬ 
tween the LE systems emergent from various symmetry¬ 
breaking phases of the Ising variables. 


A. Uniform SC and Period 2 PDW SC Phases 

We first review the uniform configuration of the Ising 
variable (and also the paramagnetic phase of the Ising 
variable), in which the SC operator will develop the same 
expectation value for all the LE system s 38 ! 49 . For the 
staggered (period 2) Ising configuration, there are two 
phases, depending on the sign of 8 J, a uniform SC state 
and a PDW state. We will solve the self-consistency 
equations for the both phases, the uniform SC state and 
a PDW state. Although the equations have the same 
form, they correspond to different phases. The case of a 
period 4 Ising configuration of the form can be 

treated in the same manner. The only difference is that 
the two phases will be an uniform SC state or a period 4 
PDW SC state. Here we will focus in the simpler period 
2 case. 


1. Uniform SC Phase 

In the uniform configuration of the Ising variable, the 
effective Josephson interaction between neighboring A 
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subsystems (the LE liquids) is 

"Hint = -Jt ^2 cos[v / 27t(6'a, j - 0 A ,i+i)\ (3.1) 

i 

To perform the interchain MFT, we consider only the 
terms involving the i-th type- A system among 'Hint- Us¬ 
ing standard interchain MF T 34 i 38 i 40 we can approximate 
eq. (13.11) by: 



with 2/i = Jt[{ cos(-\/27r0i+i )) + (cos(v/27T0i-i))]- The 
self-consistency of the MFT then requires that 

(cos(^2t T0 Aii )) = -2~. (3.3) 

Jt 

Following Refs, piolfdsl the self-consistency equation can 
be solved from the following two expressions: 


(cos(v / 27r6 , A,i)) = 

( r (| + |)F(i 

l 4 \/^ 


(l + Q 7 rr(l-d/ 2 ) , 

16 sin 7 r£ r(d/2) 

fx \ (d~ 2) 


“I) 


2 sin 




d 

M d 


(3.4) 


where M, the soliton mass in the 1 + 1-dimensional sine- 
Gordon model, is related to /i by 


F(d/ 2 ) ( 2r(g/2) 

7rr(l-d/2) ^0FT(i + |) 


d—2 


M 2 ~ d 


(3.5) 


In these equations d = 1/(2K C ) is the scaling dimension 
of the vertex operator e l '^ 0A ’ i and £ = -^iTr Using 
equations Eg. (13.41) and Eg. (13.51) . we can compute explic¬ 
itly the value of (cos(v / 27 rd A ,i)) for a given value of Jt 
and K c . This completely determines, at least at the mean 
field level, the solution of the coupled LE system s 34 ' 40 . 


2. Period 2 PD W SC Phase 


the MF solution doubles the unit cell, due to the trans¬ 
formation \J 7 2jtQ A ,’ii —> v27t(9a,2 i + 7r acting only on the 
even sites. Thus, SC order parameter oscillates in space 

AjOf) ~ (-1)- , (cos(v / 2tt6(a)), (3.7) 

corresponding to a period-2 PDW SC state. 

Before moving onto the coexistence phase in the next 
section, let us mention what is the dependence of T c with 
SJ (or Jt, depending on the Ising configuration). We 
can think of 2/i in eq. (13.21) effectively as an external held 
due to the mean held value of m :j = (cos(v / 27r(T/)) in the 
nearest neighbor systems. We can write then 



in which hj = J(rrij+i + to_/_i) and is the conven¬ 
tional kinetic term for the Luther-Emery liquid. As we 
saw above, for the uniform or staggered configuration the 
value of uij is the same in all the systems, or effectively 
the same for SJ < 0 since we can perform a transforma¬ 
tion on the even sites y/2nd A ^i — > V2nd A ^i + 7r. 

In summary, we can write just m = m,j = (cos(-\/27r(7/)} 
and h = hj = 2 Jm (where J = 8J or Jt depending on 
the case). For h —► 0 we have that self-consistency implies 

m = xsch = 2Jxscm , (3.9) 

which has the trivial solution m = 0 or a non-trivial solu¬ 
tion to J 0 if 2 Jxsc = 1 (which determines the critical 
temperature). Using that for a Lutter-Emery liquid 

XsciT) ~ T2 A 1 S /Xe , (3.10) 

we have that: 


T c ~ A s J a (3.11) 


where the exponent is a 


1 

2-1 /K c ' 


Although the re¬ 


sulting T c is small when J is small, what is important 
is that is only power-law small, instead of exponentially 
small as in the BCS case. 


In the staggered configuration of the Ising variable, the 
interaction term between the A systems is 

Hint, = ~SjJ2 COS[V^(0A,i - e A ,i+ 1)], (3.6) 

i 

which is identical to that of the uniform configuration 
case, Eq. (ED, if SJ > 0. Hence if SJ > 0, we can 
simply replace Jt by SJ to find the MF solution. This 
will give a uniform SC state. 

If SJ < 0, then we can perform a transformation on 
the even sites, \f 2 Jd A , 2 i —> VJi r$A, 2 i + 7r, effectively 
changing the sign of SJ and coming back to the first case. 
Though the form of the equation is identical to that of 
the uniform SC state, it is important to remember that 


B. Uniform SC and Period 4 PDW SC state 
coexistence phase 

Now we consider the period 8 states of the Ising vari¬ 
ables <7i = (—1) L*/ 4 J. Then the Josephson coupling also 
modulates in space with period 4, and thus we need to 
solve four coupled self-consistent equations in MFT. The 
effective MF Hamiltonian for each A system is given by^ 

Hf>l = -2/q J d 2 xcos(V2ir9 A j) (3.12) 

with 

2 Hi = [Ji(cos(V2n9 A ,i+i)) + j7 i _i(cos(v / 27r6'A,i-i))] 

(3.13) 









where Ji = Jaa~ J'aa 0 ^^^ i n which a i is i n the period 
4 structure. 

Using J T = Jaa + J'aa and 5 J = Jaa - J'aa and 
defining m, = (cos(-\/27r 6a, i), it is clear that we need to 
solve only for the four systems i = 0,1, 2,3 in this MFT 
by assuming that the MF solution does not break the 
translational symmetry i ~ i + 4 of the pattern of the 
Josephson coupling. 

Upon implementing the MFT analysis from the pre¬ 
vious section we have the following set of coupled equa¬ 
tions: 


mo = f(d) 


mi = f(d) 


m 2 = f(d) 


m 3 = f{d) 


^ m 3 SJ + TOiJr ^ d/(2 “ d) 
^ m 0 J T + m 2 J T ^ d/(2 ~ d) 
^ toiJt + m 3 J T ^j d/ ^ d) 
^ m 2 J T + d) 


(3.14) 


In the limiting case where Jt = SJ (i.e. J'aa = 0) 
Eq. (13.181) has the trivial solution x = 1. In this case all 
the SC order parameters are in phase in the case SJ > 0. 
On the other hand, for 6 J < 0, there is a shift of n every 
four lattice spacings. So, in this case, the periodicity of 
the PDW order parameter will be eight (and not four), 
although the self-consistency equations actually will take 
the same form. 

For now we will assume 5J > 0 (see section II V Cl for 
the 5J < 0 case). Then in the pattern that we consider 
here, we find x < 1 and so there is a coexistence between 
the uniform SC and the period 4 PDW order parameters. 
Let us now solve the system of equations (I3.16D - (I3.17D nu¬ 
merically for some values of the parameters. The results 
are summarize in table [T| 

We now compute T c for this case. Following the same 
steps as in the previous section we have that: 



where f(d) is a constant that only depends on the scaling 
dimension d = The explicit expression for f(d) is: 


(i + g)7rr(i - d/ 2 ) ( r( d + f)r(i 

16sin7r£ F(d/2) y 4-^7? 



x 



/7rr(l - d/2) 

V r(d/ 2 ) 


d/(2-d) 


and 

ho = SJm 0 + Jrmi 

hi = Jrmo + Jxm 2 (3.20) 

where we have used that mo = m 3 and mi = m 2 . Since 
all the A-systems are equivalent, they have the same SC 
susceptibility Then, the self-consistency equations are 


( 2 r(g/ 2 ) \ 


(3.15) 


Notice that the system of Eqs. (13.141) is non-linear. Nev¬ 
ertheless it is easy to see that mo and m 3 (mi and m 2 ) 
will take the same value (mo = m 3 and mi = m 2 ). We 
can therefore reduce Eq. (13.141) to a system of only two 
coupled equations: 


crA\( m ^J + miJT \ d,{ - 2 d) lo 

m 0 = f{d)[ ---J (3. 

' m 0 , Jt + miJ T \ £i/(2 ~ d) 


mi = f(d) 


16) 
(3.17) 


Taking the ratio of Eq. (13.161) and Eq. (13.171) we get: 

\x + 1\ 

Tfr) <3 - 18) 


where A = 8 J/ Jt- 

We can solve numerically the previous transcendental 
Eg. (13. 141) . or directly solve the system Eqs. (I3.16I) - (I3.17I) . 
Before solving the system of equations (I3.16D - (I3.17D nu¬ 
merically for some values of the parameters, let us com¬ 
ment on Eq. (13.181) . 


m 0 = Xscho, mi = xschi (3.21) 

We can write this as a system of linear equations, 

/ 1 - Xsc5j - XScJt \ fm O \_f0\ 

^ - XscJt 1 - XscJt ){mi) ^ 0 J 

which has a non-trivial solution if and only if the determi¬ 
nant of the 2x2 matrix is zero. This gives us a quadratic 
equation in for xsc- Choosing the positive solution we 
find that the critical temperature for the coexisting state 
is 


Tr = A., 


2J t (Jt-8J) 


-Jt -5J+ yjhj% - 2 J t 5J + SJ 2 J 

(3.23) 

where we recall that the exponent is given by a = 

- —Notice that, in the limit <5J —> Jt, we re- 

2 1 / K c 

cover Eq. (13. IIP Thus, as in the uniform or pure period 
2 PDW state, T c has a power law behavior in Jr and 
SJ, and it is not exponentially small as it would be in a 
weak coupling BCS type theory. 














































9 


Jt SJ d mo mi m 

niPDW 

1 

1 1/4 0.890893 0.890893 0.890893 

0 

1 

0.8 1/4 0.876601 0.889789 0.883195 

0.0065943 

1 

0.5 1/4 0.853007 0.887947 0.870477 0.0174703 

1 

0 1/4 0.806035 0.884205 0.845120 

0.0390853 


TABLE I. Numerical solution for the system of equations 
(13.1711 for different values of the parameters Jt, dj and 
d — 1/2K C . We also define m = (mi + mo)/2 and mpDW - 
(mi-mo)/2, which correspond to the uniform and PDW part 
of the SC order parameter. 

IV. FERMIONIC QUASIPARTICLES OF THE 
SUPERCONDUCTING STATES 

So far, we have solved the coupled LE systems in the 
limit \Jaa,]\ > \Jab\ and \Jaa,j\ » \J'ab\ in Eq.(EH) 
so that the couplings of the LE systems to eLL systems 
can be taken as the perturbation. In this limit, we have 
ignored the type- B eLL systems and shown that the var¬ 
ious SC states can emerge. Now we include the eLL sys¬ 
tems back and investigate the nature of the full emergent 
SC state by looking at the SC proximity effect. First of 
all, we note that the eLL systems themselves will flow to 
the 2D Fermi liquid fixed point (at low enough temper¬ 
atures) under the effect of the hopping amplitude tsB- 
This is the most relevant coupling in Ea. (12.9l) . The re¬ 
sult, for t bb small enough, is an anisotropic Fermi liquid 
with an open Fermi surface, shown as the dashed curves 
in Fig (21 

Having solved the largest energy scales in Eg. (12.91) . set 
by tBB and Jaa, we now include the effect of the pair 
tunneling processes mixing the systems A with the sys¬ 
tems B 1 presented in Ea. (l2.9D . and parametrized by the 
coupling constants Jab and J' AB , respectively. We will 
study the effects of the SC states on the A systems on 
the B systems by treating the pair-tunneling terms to the 
lowest non-trivial order in perturbation theory in these 
coupling constants. Hence, we are assuming that the in¬ 
teraction with the type-H eLL systems does not back re¬ 
act to considerably change the MFT value of the SC gap 
in the LE systems. As in Ref. j3f|, under the proximity 
effect mechanism the B systems become superconduct¬ 
ing and provide the quasiparticles for the combined A-B 
system. 

Since we are interested in the effect of the SC order pa¬ 
rameters on the electronic spectrum, we replace the pair 
density A A,j(x) of the type-A LE systems in Eg. (12.91) by 
its MF value (A A,j) determined by the interchain MFT 
discussed in the previous Section m In this approxima¬ 
tion, we find that Eg. (12.91) reduces to 

H' -A ^ / dx { ~ t BB'^2[lpB Jta i>B,j+l,a + h.C.] 

j a 

— Jab[Al b ^{/S.A,j) + Agj(Ayjqi) + h.C.] 

++ V'Bj-l.tV’Byq) +h.c.]|, 

(4.1) 


which is simply a theory of a Fermi surface coupled to the 
SC via a proximity coupling. Since Eq. m is quadratic 
in the electron fields, we can readily diagonalize the effec¬ 
tive Hamiltonian, and obtain the quasiparticle spectrum 
for the different SC states found in Section m 


A. Uniform SC phase and pure PDW phase 

As we saw in section IIII A 21 for the staggered (period 
2) configurations of the Ising variables is possible to have 
either a pure uniform SC state or a pure PDW state. The 
case of the uniform SC was studied by Granath, et al 
who showed that, depending on the values of Jab and 
Jab i it is possible to have either a d-wave SC state with a 
fully gapped spectrum of quasiparticles or a conventional 
d-wave SC state one with a nodal quasiparticle spectrum. 
We refer the reader to their paper for further details^ 

On the other hand, for the pure PDW state, even 
though the MF equation for the SC gap has the same 
form as for the uniform SC gap, the quasiparticle spec¬ 
trum is quite different. We will study this spectrum in 
detail here. Let us start by defining the period 2 PDW or¬ 
der parameter, i.e. with ordering wave vector Q = (0,7r), 

A f = A Q e™ j (4.2) 

where Aq is given by the spin gap and the interchain 
MFT value for (cos \Z2tt9) which is given in Section UlI Al 
for the period 2 configuration of the Ising variables. No¬ 
tice that for a period 2 state A q = A_q, since for a 
period 2 state Q and Q differ by a reciprocal lattice 
vector. 

To find the quasiparticle spectrum we first write down 
the Hamiltonian of Eq. (E3D in momentum space. Defin¬ 
ing the Nambu basis (here we dropped the B label in the 
electronic operators, since it is understood that we are 
referring to the eLL systems) as: 

*k = ^k+(0,7r)t> ^- k P V>—k—( 0 ,ttR) (4-3) 

we can write the Bogoliubov de-Gennes (BdG) Hamilto¬ 
nian as 

H = Y, ^k ^k V'k (4.4) 

k 


where the one-particle Hamiltonian H k is given by 
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’e(k x ) - tsB cos (k y ) 0 

0 e(k x ) +t B BCOs(k y ) 

0 2iJ' AB /S. Q sin (k y ) 

\ — 2iJ' AB &.Q sm(ky) 0 


From this one-particle Hamiltonian we find the quasipar¬ 
ticle spectrum 


E( k) = ±t B B cos (k y ) ± \Je 2 {k x ) + 4J^ 2 b |Aq| 2 sin 2 (fc !/ ) 

(4.6) 

In Fig. [2] we plot the Fermi surface of the Bogoliubov 


J V 

C 0 

A A 


-n -- -- 0 - - 3n tt 
4 2 4 4 2 4 

k x 

FIG. 2. (Color online) On the left, Fermi surface for the pure 
period 2 PDW state. The dashed (blue) line corresponds to 
the original FS in the absence of superconductivity. The solid 
(red) line corresponds to the new FS after the superconduct¬ 
ing proximity state is established. On the right, the spectral 
function T(k, 0) corresponding to the pockets on the left. We 
used J' ab Aq = 0.12t, tsB = 0.7 1, e(k x ) = —tcosk x and 

6 = 10~H 


\T 

V/ 

A 

A 

\) 

\/ 

V 

A 

A 

A 


4 2 4 4 2 4 

K 


quasiparticles of this period 2 PDW state for some val¬ 
ues of the parameters. In contrast to the pure uniform 
SC state, whose spectrum can be either nodal or fully 
gapped, we find that this PDW state (Aq ^ 0 in Eq. 
(14.2D ) has pockets of Bogoliubov quasiparticles, as it is 
also found in the weak coupling theories ^ 14 ' 17 ! 44 ; 50 ! 51 
The size of the pockets depends on the strength of the 
SC gap. In addition, we compute the spectral function 
given by (see for instance Ref. [52|): 

H(k, w) = —-Im[Gn(k, w)] (4.7) 

7r 


where 


G( k,w) =-^ 

w + id — Hk 


(4.8) 


is the retarded Green function and <5 = 0 + . The spectral 
function H(k, ui = 0) for this pure period 2 PDW state is 
shown in Fig. [2j In Fig. [3]we plot the dispersion relation 
of the Bogoliubov excitations for several values of k y . 


0 

—2iCf' A B^*Qsin{ky) 
-e{k x ) + t B B cos (k y ) 

0 


^J'ab^q sinfe) \ 
0 
0 

-e(k x ) - t B B cos (k y )J 


(4.5) 



FIG. 3. (Color online) a-d: Dispersion relation of the Bo¬ 
goliubov quasiparticles for k y = 0, j, respectively, for 

the period 2 PDW state. Here we used J' ab IS.q = 0.12t, 
t B B = 0.71 and e(k x ) = —tcosk x 


B. Coexistence Phase of a Period 4 PDW and a 
uniform SC: the Striped Superconductor 

We start by writing the SC order parameter, which 
includes both the uniform SC and the PDW order pa¬ 
rameters as an expansion of the form 

Af = A 0 + v^Aq cos ^ (4.9) 

where the expectation values of the order parameters Ao 
and Aq, where Q = (0, -|) is the ordering wave vector, 
are set jointly by the spin gap of the LE systems and 
by the interchain MFT value for (cos \/2iiff) found in the 
previous section for the period 4 state of the Ising degrees 
of freedom. 

We now write down the Hamiltonian in momentum 
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FIG. 4. (Color online) Quasiparticle spectra with nodal 
points in the coexistence phase. The dashed (blue) line corre¬ 
sponds to the original FS in the absence of superconductivity. 
The red points correspond to the position of the nodes in the 
absence of the PDW state (A q = 0). The green points cor¬ 
respond to the position of the nodes in the presence of the 
PDW state with A q = 0.2. We have chosen the parame¬ 
ters J' AB = 0.5f, Jab = 0.2 1, tBB = 0.7 1, Ao = 0.2 and 
e(fca) = —tcosk x 
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space following the notation of Ref. [50}. We define the 
Nambu spinor as: 

= (V’k-p ^k+qt> ■ • ■ ’ V’-kf'. ^-(k+q)4., • • •) (4.10) 

where q is the ordering wavevector. In our case q = 
(0, 7 t/ 2 ) and k is taking values over the reduced Brillouin 
zone (RBZ) associated with the ordered state, which in 
this case is k x £ [— 7 r, 7 r) and k y £ [— 7 r/ 4 , 7 t/ 4 ). In this 
basis the Hamiltonian is given by: 

H = J2 Wc^k, (4.11) 

ke/tsz 

where the BdG Hamiltonian in the Nambu basis of 
Eg. (14.101) is given by: 

Md-l) (412) 

where Ak = diag(e(k), e(k + q),...) is a diagonal ma¬ 
trix, and the square matrix Ck contains the SC order 
parameters. Since the ordering vector is 7 t /2 along the 
k y direction, our matrix Ck is given by a 4 x 4 matrix 
with the form: 

/ / 0 (k) A(k) f 2 ( k) / 3 (k) \ 

c = A*(k) /o(k + q) /i(k + q) / 2 (k + q) 

k / 2 *(k) A*(k + q) /o(k + 2q) /i(k + 2q) 

V/|(k) fU k + q) /i (k + 2 q) / 0 (k + 3q)/ 

(4.13) 


where /o corresponds to uniform pairing and f \, / 2 , / 3 to 
the finite momentum pairing. The explicit expressions 
are the following: 

/o(k) = 2A 0 (J ab - Jab cosk y ) 

/i(k) = - iA q (Jab - V2, Jab cos (k y + q y / 2)) 

M k) = 0 

/ 3 (k) = iAq(J A b - V2Jab cos{k y - q y / 2 )) (4.14) 

where we recall that q = (0, 7 t/2 ), so q y = 7r/2. 

First of all, due to the periodicity of the PDW SC 
state, it is necessary to the fold original FS. Let us first 
analyze the case of the pure uniform SC state. In this 
case (Aq = 0) the spectrum can be easily calculated 
from the Hamiltonian given in eq. (14.121) : 

El ± = ( e(k x ) ± t B B cos (k y )) 2 + A 1 (Jab ± J'ab cos (k y )) 2 

El ± = ( e{k x ) ± t B B sin(fc !/ )) 2 + A 1(Jab ± Jab sin ( fc y )) 2 

(4.15) 

We can see that this SC state will have a quasiparticle 
spectrum with nodes if \Jab\ < \J'ab\- Now, even in the 
coexistence phase, where both Aq A 0 and Ao A 0, the 
quasiparticle spectrum may still can have nodes. For the 
pure uniform SC state, the position of the nodes depends 
on the values Jab I J'ab an< 4 ^bb- In the coexistence 
phase the position of the nodes will depends on Aq as 
well (see Fig. a. As in the case of pure period 2 PDW 
state, we show in Fig. [5] the dispersion relation of the 
quasiparticles for several values of k v . 


C. Period 8 PDW state 


Above we focused on the coexistence phase for the pe¬ 
riod 4 case. This was the case when 5J > 0. However, if 
SJ < 0 (i.e. for Jaa < J'aa) the case is different and we 
find a PDW state. There is a shift of n every four lattice 
sizes, so in this case the periodicity of the PDW order 
parameter is actually eight (not four!). Nevertheless, the 
self-consistency equations will have the same form: 


m 0 = f(d) 


mi = f(d) 


^ wqI^I + d) 

^ttioJt + rriiJr^j ^ ^ 


(4.16) 


The pattern of the SC order parameter is now that of a 
pure period 8 PDW SC state: 


A = (Ai, A 2 , A 2 , Ai, — Ai, — A 2 , — A 2 , — Ai, Ai,...) 


We can write the previous pattern using the following SC 
order parameter: 


A 3 4 = A sin 



+ A sin 



37T 

~ 8 ~ 


(4-17) 
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FIG. 5. (Color online) Dispersion relations of the quasi¬ 
particles in the coexistence phase shown for several values 
of k y . Notice that the dispersion relation is only gapless 
for k y « 0.457, which corresponds to the position of the 
nodal point for the same set of parameters used in Fig. []] 
J'ab = 0.5t, Jab = 0 .2t, tBB = 0.71 , Ao = Aq = 0.2 and 
e(k x ) = —tcosk x 



3JT 7T () H H 37T jjr 

~~4~ ~2 4 4 2 T” 


k x 

(a) 



8 _7T 3JT _7T Q 7T 7T jj 

~ 4 ~2 4 4 2 T 


k x 

(b) 



8 _jC 3n tt 0 — — 377 jj 

~~A ~2 4 4 2 T 


k x 

(c) 

FIG. 6 . (Color online) FS for the period 8 PDW state. The 
dashed (blue) line corresponds to the original FS in the ab¬ 
sence of superconductivity. The solid (red) line corresponds 
to the new FS (pockets) after the introduction of supercon¬ 
ductivity (Ai 7 ^ 0 and A 2 7 ^ 0). In (a) Ai = A 2 = 0.05, 
in (b) Ai = 0.08 and A 2 = 0.1 and in (c) Ai = 0.25 and 
A 2 = 0.3 In all figures J' AB = 0.6 1, Jab = 0.4t, tBB = 0.7 1, 
Ao = 0 and e(k x ) = —tcosk x 


where we have defined: 

A = Ai sin Q) + A 2 cos 

A = Ai cos - A 2 sin (0 (4-18) 

where Ai and A 2 are given by the spin gap and the 
interchain MFT value for (cos\/27t 0) in eq. (14.161) . 

Since we are dealing with a period 8 SC state, the re¬ 
duced Brillouin Zone is now case is k x £ [—w, n) and 
k y £ [— 7 t/ 8 , 7 t/ 8 ) and q = (0, 7 t/ 4 ). The difference be¬ 
tween the period 4 and the period 8 is that the definition 
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of the Ck matrix is different, since is now an 8 x 8 matrix. 



f /o(k) 

/i(k) • 

• A(k) \ 

Ck = 

A*(k) 

/o(k + q ) ■ 

• /e(k + q ) 


Ur( k ) 


/o(k + 7 q ) J 


(4.19) 


where the /i(k)’s are given by the following expressions: 

/o(k) =/ 2 (k) = / 4 (k)=/ 6 (k) = 0 

/r(k) =iA Q J AB (e-™' 8 + e"^ 8 ) 

-J' A B e ~ m/A cosiky + g y / 2 )) 

/a(k) =*A Qj AB (e- 3 ^ + e"' 9 ^ 8 ) 

-J' A B e ~^ /A cos (ky + 3q y /2)j 

M k) = - *A Q JA B (e 3i7r/8 + e i9 "/ 8 ) 

-J A B e ^ /A cos(ky - 3q y /2)) 

/ 7 (k) = - «A Q J AB (e i7r/8 + e 13 ^ 8 ) 

-J AB e ™ /4 cos(fc y - g y / 2 )) 


where we recall that q = (0,7r/4), so q y = 7 t/4 . Having 
Ck we can write down our BdG Hamiltonian as in eq. 
(14.121) . In Fig. ([ 6 | we show the FS for some values of Ai 
and A 2 . As in the pure period 2 PDW state, we see the 
formation of pockets due to the folding of the FS. 


V. OTHER PHASES 

For completeness we summarize the other possible 
phases occurring in the system. Following closely 
Granath et al — we treat the interactions appearing in 
eq. (12.91) perturbatively around the so called decoupled 
fixed point. At this fixed point (FP) the systems are 
completely decoupled, and each one of the systems corre¬ 
sponds to a ID system that can be solved using bosoniza- 
tion. Around the decoupled FP a perturbation with cou¬ 
pling constant g is relevant (irrelevant) if its scaling di¬ 
mension d g < 2 ( d g > 2). The scaling dimensions for the 
operators appearing in Eq. (12.91) are given in the work 
of Granath et al.J^ The phases found by Granath et al. 
are: 

1. Typically, the couplings between the eLL and LE 
systems are irrelevant or less relevant than the cou¬ 
pling between AA and BB systems separately. In 
this case the RG flows to the point where all the 
AB couplings go to zero. At this FP the system 
is made of two (independent) interpenetrating sys¬ 
tems, A and B. 


2. The J AA ( J ' AA ) term is relevant for > 1/2. 

In this case the A systems develop long-range order 
and a full spin gap. Since the BB electron tunnel¬ 
ing operator has lower scaling dimension than the 
BB spin exchange interaction, in the absence of 
a charge gap in the B subsystem, most probably 
the B subsystem is in a anisotropic Fermi liquid 
phase. However, this two fluid FP is unstable due 
to the proximity effect. Depending on the parame¬ 
ters in the Hamiltonian of Eq. (12.91) the quasiparticle 
spectrum can be gapless (present nodes or pock¬ 
ets in the pure PDW state) or fully gapped. This 
means that we can have several possible stable SC 
phases, a SC state with Fermi pockets, a nodal SC 
state, or a fully gapped SC state. These were the 
phases studied in the previous sections using inter¬ 
chain MFT and coupling the eLL systems to the 
LE systems. 

/ D\ 

3. If Ac > 0, the B subsystem can develop a an¬ 
tiferromagnetic phase. At this FP will be a coex¬ 
istence between superconductivity (in the A sub¬ 
system) and antiferromagnetism (in the B subsys¬ 
tem). This FP is stable, due to the spin gap in the 
SC (A) and the charge gap in the antiferromagnet 
(. B ). The quasiparticle spectrum is therefore fully 
gapped as is also found in BCS-type theories.— 


VI. CONCLUDING REMARKS 

We have investigated a model of an array of two in¬ 
equivalent systems in the quasi-one dimensional limit. In 
this limit we have treated the interactions between the 
different systems in the array exactly using bosonization 
methods and interchain mean field theory. The phases 
that we found are either a uniform d-wave superconduc¬ 
tor, a striped superconductor (in which the uniform SC 
and the PDW SC state coexist), and a PDW state. To 
simplify the analysis we only looked at the case in which 
the modulation of the SC state is commensurate. 

The resulting critical temperatures are, as expected, 
upper bounds on the actual phy sical critical tempera¬ 
tures. As emphasized in Refs. [40] and [53j, the analytic 
dependence of these mean field T c ’s on the coupling con¬ 
stants obeys the exact power-law scaling behavior pre¬ 
dicted by a renormalization group analysis of the di¬ 
mensional crossover from the ID regime to the full (but 
anisotropic) 2D phases, albeit with an overestimate of 
the prefactor. 

On the other hand, the actual critical temperatures are 
significantly suppressed from the values quoted here due 
to the the two-dimensional nature of the array. Hence we 
expect the ground states that we found here to undergo 
a sequence of thermodynamic phase transitions leading 
to a complex phase diagram of the type discussed by 
Berg et al — (and by Agterberg and Tsunetsugu4) It 
is well known from classical critical phenomena of 2D 
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commensurate systems that states of the type we discuss 
here may become incommensurate at finite temperatures 
due to thermal fluctuations if the period of the ordered 
state is longer than a critical value (typically equal to 
four), see, e.g. Ref. (H3| . 

We have shown that a high energy scales (of the order 
of the spin gap), we can first determine the SC phases 
of one set of systems (in our notation, the Luther-Emery 
liquid systems A). At these energy scales we showed that 
it is possible to have, in addition to a uniform SC phase, 
a pure PDW state and a coexistence phase of a uniform 
and a PDW state. Having determined the SC in the 
LE systems, we proceeded to incorporate the electronic 
Luttinger liquid systems perturbatively. We found that 
the quasiparticle spectrum arising from the eLL systems 
can present Fermi pockets if the SC state is a pure PDW 
state. In the case of coexistence uniform SC and PDW 
state or pure uniform SC (i.e. a striped superconductor) 
the quasiparticle spectrum can have nodes or be fully 
gapped depending on the value of the coupling in the 
model. We should stress, as it was done recently in Ref. 


GU , that in this quasi-ID approach the superconducting 
state evolves from a local high energy scale, the spin gap, 
which hence has magnetic origin. For temperature T 
higher that the spin gap, the system is a quasi ID system 
which does not have quasiparticles in the spectrum up to 
a scale, determined by an electron tunneling scale, to 
a crossover to a Fermi liquid type system. Hence, at 
least qualitatively, systems of this type behave as ‘high 
T c superconductors.’ 
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